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Abstract 

The peculiar features of the mixing in the neutral pseudoscalar mesons K° — K° can be 
introduced by the analogy to the optical polarization. The time-reversed not-invariant 
processes and the related phenomenon of CP-nonconservation can be then joined 
to the dissipative effects which yield a not vanishing imaginary part in the relevant 
propagation of electromagnetic radiation in a medium. Thus, the propagation of the 
two transverse polarization states can reproduce the peculiar asymmetries which are 
so common in the realm of high energy particle physics. 



1 Introduction 



Mechanical models that provide a viable description of the non-unitary time evolution 
can be related to the problem of two coupled degenerate oscillators with dissipation Q. 
The physical intuition underlying any proposed model stems back from the fact that 
any time-reversed non invariant process can yield a non vanishing imaginary part for the 
relevant Hamiltonian. Although, the Hamiltonian of a complete sensible quantum system 
is expected to be a Hermitian operator, under suitable conditions we may recover the 
time evolution of a subsystem according to an effective non hermitian Hamiltonian like 
in the case of metastable states HJ. A celebrated example where this description has 
proved extremely useful is the two-states kaon complex |J. In this case, the single pole 
approximation of the Weisskopf-Wigner method jEj can be applied to derive the eigenstates 
Kg and Kl of a 2 x 2 effective Hamiltonian TL [^J. Although such puzzle system requires 
the formalism of density matrix, and the notion of the rest frame for an oscillating unstable 
composite system appears difficult to be implemented, nevertheless there has been always 
a steady activity towards understanding these processes at the level of wave- function Q. 
The time evolution of the flavour states K° and K° can be written as 



lh dt{\K°(t)))- n {\K*(t))) 



(H n H l2 \f\K^_(t))\ 
\H 21 H 22 )\\K\t))) 



(1) 



where the effective 2x2 matrix Hamiltonian TL is linear but not necessarily Hermitian and 
it can be uniquely rewritten in terms of two Hermitian matrices M and V: TL = M — iT/2. 
M is called the mass matrix and T the decay matrix. Their explicit expressions can 
be derived by the weak scattering theory responsible of the decay Q. The dynamical 
behaviour of this system can be described in the Schrodinger picture by means of a time 
evolution operator O: 

which can be written by means of the well known exponential solution 
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0(t) 

whose matrix elements are given by 



exp 
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where as usual the latin indices denote K°, K°. In order to evaluate this matrix expo- 
nential form, we can use the Sylvester's formula of the matrix spectral decomposition for 
a function f(A) of a generic n x n square matrix A. It is easy to express this Sylvester's 
formula for the special case of 2 x 2 matrix with two distinct complex eigenvalues a% and 
a 2 
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d 



f{A) - -Q!^ 



1 + 



df(ai) 
da i 



A = f{a 1 )l + f{a 1 ){A-a 1 I) 



(5) 



(6) 



2 



This formula becomes rewarding in connection with the expansion of the exponential 
representation of the evolution operator. In fact, in this case, we get 



0(t) =exp ( ~Ht] = - 

n J ai- 
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It is convenient to represent the matrix TL in terms of the Pauli spin matrices 

ft = h^a^ = h l + h • er 

where we let 



(7) 
(8) 
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being trH = a\ + 02 = 2/io, detH = aia2 = (ho) 2 — |h| 2 and then 02 — ot-i = 2|h|. It is 
easy then proven that if is a normal matrix with a complete set of orthogonal eigenvectors 
if and only if h x h* = 0. Rewriting the solution in terms of the h vector, we have 
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where h = h/|h|. Except for the overall phase exp (—ihot), this is just a rotation 7£[u;t] = 
exp [—i(uj ■ ct/2)] of the initial state |-0(O)) : 



|^(t)) = e-^ftMIV^O)) • 



(11) 



In particle physics, K° and K° are expected to be distinct particles from the point of view 
of the strong interactions, they could transform into each other through the action of the 
weak interactions. In fact, the system K° — K° results degenerate due to a coupling to 
common final states (direct CP-violation) (K° <-> -kit <-> K°) or by means of a mixing 
(K° <-> K°) (indirect CP violation). Even if this is an effect of second order in the 
weak interactions, the transition from the K° to K° becomes important just because this 
interference in evaluating decay widths. The mass eigenstates are linear combinations of 
the flavour states |K°) and |K°) namely \Kg > and >, which have definite masses 
and lifetimes. A convenient way to discuss the K° — K° mixing problem is then to express 
Ks and Kl as eigenstates of the 2x2 effective Hamiltonian. It can be shown that if CPT 
holds then the restriction Tin = ft-22 and hence Mn = M22, ^11 = ^22 must be adopted. 
Furthermore if CP invariance holds too, then besides Hu = H22 we get also H12 = H21 
and consequently Ti2 = I^i = T2i*, M12 = M 2 i = M21* so that r^, Ma will result all 
real numbers. The time evolution of the mass eigenstates 



\K L (t))) V[t) 



\Ks)\ 
\Kl)J 



(12) 



is governed by the matrix elements 

V a p = (K a \ exp 

where greek letters denote K$, Kl- The evolution matrices U and V are then related by 
the following similarity transformation 

U = UVR- 1 (14) 

where the complex scaling 1Z in any CPT invariant theory can be found by means of the 
the eigenvalues of the effective Hamiltonian TC 



As = H u - VH 12 H 21 = M n - -T n -Q = ms --T s (15) 
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X L = H n + VH 12 H 21 = M n - % -Tu + Q = m L - % -T L (16) 



where 



Q = Vh 12 h 21 = J(m 12 - -r 12 )(M 12 * - -r 12 *). (17) 



These real (ms^) and imaginary (Ts,l) components will define the masses and the decay 
width of the TL eigenstates Ks e Kl . These short- and long-lived particles result then a 
linear combination of K° and K°: 

where usually TZ l is parameterized according to the following relation 

\P q J y/T+W Via/ 
being \p\ 2 + \q\ 2 = 1. With a proper phase choice, the mixing parameter 



Q H 21 

a 



p \ h 12 \ m 12 - |r 12 



M* —IV* 
12 2 12 (20) 



manifests its presence in the non orthogonality between the oblique independent states 
Kl and Ks- 

(K S \K L )= 1 —^ . (21) 
1 + |a| 

The experimental evidence in 1964 that both the short-lived Ks and long-lived Kl 
states decayed to irir upset this tidy picture. It means that the states of definite mass 
and lifetime are never more states with a definite CP character. With the conventional 
choice of phase, the CP eigenstates K\ and K 2 enter into play so that the mass eigenstates 
can be parameterized by means of an impurity complex parameter e which encodes the 
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indirect mixing effect of CP violation in the neutral kaon system. The general expression 
of this impurity parameter 

e*? /ImMi2 i ImFi2 
6 = 271 VReM 12 ~ 2ReM 12 

where £o represents the additional manifestations of the CP violation due to the effective 
isospin decay amplitudes. Neglecting these decay effects, to derive the amount of this 
complex parameter, it appears more convenient to write 

q 1 — e 

a = - = - 

p 1 + e 

and then 

_p — q _ 1 - a _ \/Hi2 - VH21 
p + q l + a VH12 + VH 2 i 
_ 2iImM 12 + Imr 12 

~ (2ReMi2 - iRel^) + (Am - §Ar) (Am - §AT) 

Thus it results evident that the CP-violation parameter e arises from a relative imaginary 
part between the off-diagonal elements M\2 and Ti2 i.e. if arg{M*^T\2) 7^ 0. Although 
the discovery of CP- violation || indicated that the kaon system is somewhat more com- 
plex than the typical oscillating two-state problem and involves considerably more subtle 
complications; nevertheless there are many attempts to illustrate this system in classical 
physics 0. In this context, the origin of the space-time discrete symmetries and their 
violation represent one of the most controversial topic in the realm not only of physics. 
For a long time, the fact that Maxwell equations were invariant under space-inversion 
or parity (P) and time-reversal (T) bolstered the idea that all the laws of physics are 
invariant under those discrete operations. It was easily seen that electromagnetic equa- 
tions possess another discrete symmetry since they are unaffected by a charge conjugation 
(C) operation which reverses the sign of all the charges and converts a particle into its 
antiparticle. However, since 1964 ||, we know that CP is also violated (although to a 
much lesser extent than the amount of parity violation ||) in weak interactions among 
fundamental particles. On the other hand, the origin of CP-violation is still not explained 
since CP- violating tiny effects are known to be smaller than the usual weak interaction 
strength by about three orders of magnitude and it is not excluded that CP-violation 
could be an indication of some effects of new physics at a much higher energy scale. The 
only almost overwhelming theoretical prejudice comes against CPT violation. There are 
very strong reasons || to believe that fundamental interactions can never violate CPT 
invariance. Thus, except some loosing theoretical models, the validity of CPT is assumed 
and consequently the T violation is supposed to be reflected immediately in a CP coun- 
terpart. However, it should be borne in mind that observation of a T odd asymmetry 
or correlation is not necessarily an indication of CP (or even T) violation. The reason 
for this is the anti- unitary nature of the time reversal operator in quantum mechanics. 
As a consequence of this, a T operation not only reverses spin and three- momenta of all 
particles, but also interchanges initial and final states. Put differently, this means that 
final-state interactions can mimic T violation, but not genuine CP violation. Presently, 
genuine experimental evidence of CP-violation comes from the mixing and decays of the 
neutral kaon system. In fact, the intrinsic dissipative nature of this unstable system and 
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) = (ImM 12 + 2£ ReM 12 ) (22) 
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ImM 12 + Im(r 12 /2) ! 1 ' ' 



its decay sector, faced with the problem of the complex eigenvalue of the Hamiltonian and 
therefore with the extension of the Hilbert space which is a common tool to deal with 
open systems far from equilibrium. 

In general, we can resume that there are three complementary ways to describe the 
evolution of the complex neutral kaon system: 



1) In terms of the mass eigenstates Kl,s, which do not posses definite strangeness 

\K s (t) >=\K S (0) > e~ iXst 
\K L (t) >=\K L (0) > e~ iXLt 



t\ s = im s + — 



i\l = vniL + — • 



(25) 



2) In terms of the flavour eigenstates, whose time evolution is more complex 



\K\t) >=/ + (t)|K°(0)> + -/_(t)|K°(0)> 

a 



with 



K"(t) >=a/_(i)|K u (0)) + / + (t)|K u (0)) 



(26) 



(27) 



3) In terms of the CP-eigenstates K\ and K2 

"--0\ , \T70\ 



|^>=-^ ( |K U ) + |K U )) 

\K° >=^(|K°) - |K )) 

which we let us express the mass eigenstates as 

1 

1 



CP\Kl >= + \K^ > 
CP\K$>=-\K°> , 



\K° S >- 



(28) 



(p + 5)|^i >+(p-g)|^2°>) 
The three bases {Ks,Kl}, {K°,K°} and {K±, K2}, are completely equivalent. 



(29) 



2 The Dissipative Effects in Coupled Oscillations. 



In classical physics, an analogue of the two-state mixing problem which leads to a non-zero 
value of e seems to strain its purport more than it is necessary, since the equations of motion 
in classical mechanics are time-reversal invariant. The main features of irreversibility 
enter only considering the effects of dissipation. Anyway, these results seem to reflect the 
well known requirements of additional complementarity relations, which occur at classical 
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level, to make the equations of motion of dissipative systems derivable from a variational 
principle |J. We can recover the theory for the non conservative systems in classical 
physics in terms of the standard formalism of analytical mechanics. In the case of n 
independent generalized coordinates q k , the equations of motions are given by 

d dT dT „ 

dtdq s ~~St = j <30> 
where the generalized momenta are of the form 

Qj=f2A j q i + B ij ^. (31) 

i=i 

In the special case of conservative system, a configuration of stable equilibrium is specified 
by the following Lagrangian equations 

d dT dU dC d dC „ ^ 

dtdq k dq k dq k dt dq k 

As usual, the kinetic tensor 

T =2^1 T jkWk (33) 

j,k 

U = ^^2^jkQjqk (34) 

j,k 

are supposed bilinear without any loss of generality. Of course, the exact form of and 
(Uij) specifies the way in which the motions of the various coordinates are coupled. The 
problem of determining the normal solutions reduces to that of simultaneous diagonaliza- 
tion of these two real symmetric matrices associated with each of the kinetic and potential 
energies. The problem consists practically in finding a coordinate transformation which 
simultaneously diagonalizes both and (Uij). This is equivalent to the following set of 
second order homogeneous differential equations 

J2(Uj k qj +Tj k qj) = . (35) 
3 

This underlying mathematical problem can be formulated elegantly in matrix language 
since it can be compactly associated with the following Lagrangian 



and the potential function 



C = T-U 



1 



q*Tq - q^q . (36) 



Then the problem to determine the normal modes reduces to an eigenvalue problem involv- 
ing the simultaneous diagonalization of two real symmetric matrices, classically associated 
with the kinetic and potential energies. The recipe for carrying out this diagonalization 
is straightforward. The kinetic and potential energy matrices are simultaneously diago- 
nalized by means of a sequence of three transformations from the original to the normal 
coordinates £ 

q = R % q ' = R % S q" = R X S R 2 £ . (37) 
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In this new normal coordinates the equations of motion are uncoupled. The net trans- 
formation matrix G = R\ S R2 is not orthogonal because S is not orthogonal. This 
matrix G represents the most general (homogeneous) linear transformation which includes 
two rotation matrices R±, R2 and a contraction S. We have to sacrifice the property of 
orthogonality in order to construct a matrix that can simultaneously diagonalize both T 
and IA. The effect of this transformation is to diagonalize the kinetic and potential energy 
matrices: 

diag T = G* T G 

(38) 

diag U =G l UG . 

To be more explicit, in the case of normal solution of the form q- = aje^-V we obtain 
n linear homogeneous algebraic equations: 

^2(Uijqj - u) 2 %j)aj = , (39) 
i 

where the common phase factor e^' - ^ has been cancelled. A non trivial solution corre- 
sponds to the secular equation 

det(Uij - u?%j) = . (40) 

In general, this n-order spectral equation in the eigenfrequencies to 2 cannot be solved sim- 
ply. Nevertheless, in the case of dissipative systems the coupled differential equations that 
describe the damped motion can be obtained only by means of the use of the generalized 
momenta in Eq. (pi]). In general we must consider three matrices %j, A%j and Bij. A 
normal mode solution of the kind qj = Qje~ lut can be derived from the secular equation 

det(r ij w 2 + Aij + iuoBij) = . (41) 

Using the modern Dirac notation, although we are discussing classical quantities, its phys- 
ical solution derives from the characteristic equation 

H\q> =u 2 \q> , (42) 

being 

H = itoB - A . (43) 
Thus the equations of the motion can be rewritten as 

-%^>-^- <«> 

This set of equations resembles the case of a system of coupled damped harmonic oscilla- 
tors [10]. The dissipative effects are parameterized by the presence of small off-diagonal 



terms in TL, as it was carefully illustrated in a recent publication |hj] in the classical frame- 
work of mechanics and by means of electromagnetic coupled circuits. The time-reversed 
not-invariant processes can be induced by dissipative effects which yield a not vanishing 
imaginary part for the relevant Hamiltonian. 
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3 Electromagnetic Constitutive Relations. 



The description of the fundamental nature of the electromagnetic interactions induced 



in a medium is obscured by several mathematical problems [11]. These difficulties are 
simply an expression of our ignorance about the dynamics of all processes going on in the 
medium or our inability to quantify them in detail. From a purely macroscopic point of 
view, the crucial question is to append certain constitutive relations in dependence with an 
effective parameterization of the induced charges and current densities and feeding them 
into the Maxwell's equations. The presence of the induced sources is taken into account 
by introducing a few empirical constants in the equations. The number of such constants 
necessary depends on the nature of the medium. For the simpler linear isotropic mate- 
rials, it is well known that one needs two scalar quantities: the dielectric constant e and 
the magnetic permeability fx, to describe all phenomena of macroscopic electrodynamics, 
according to the relations 

D = e E + P = eE , H = M = . (45) 

Once the empirical constants e and \x are known, one can then derive the electric E and 
magnetic B vector fields. However, for complex systems one needs in general more than 
these two constants to describe electrodynamical phenomena in them. In many instances 
this inadequacy may be addressed by postulating a more general set of constitutive rela- 



tions to replace those in Eq. (J45J): 

D = eE + /3B, (46) 
H = ^ 1 B + 7 E. (47) 

We should comment that the restriction of the number of electromagnetic constants de- 
pends on the nature of the medium in question, since in the most general case they are 
rank-2 tensors [jll|1 . While the approach of postulating these general constitutive relations 
are well known [11], until now, it is not deeply discussed in connection with the symmetry 



requirements which could be of great help in analyzing the dynamics of propagation. The 
symmetry requirements under spatial translation or under rotations are well-documented 
experimentally, and no physicist has serious doubts about their validity as a consequence 
of the conservation of momentum and angular momentum. Invariance under spatial in- 
version and time reflection, on the other hand, seems to hold only approximately, for some 
classes of phenomena. The point is that spatial inversion symmetry, which is associated 
with parity (P) invariance, is preserved by strong and electromagnetic interactions but vi- 
olated by the weak force || . Since the latter interaction is at least six orders of magnitude 
weaker than its strong and electromagnetic counterpart, parity violation can generally be 
neglected. Similar comments apply to the case of time reversal (T) symmetry, which is 
found [||] to be violated only by an effect even more feeble than the weak force. Maxwell's 
equations have another discrete symmetry since they are unaffected by the operation under 
which p, j, E and B all change sign. This operation, called charge conjugation (C) because 
the sign of all charges are reversed, plays an important role in quantum field theory, where 
it converts a particle into its associated antiparticle. The transformation properties of the 
electromagnetic fields, potentials and currents, with respect to the discrete symmetries C, 
P and T, and some important combinations like CP and CPT (Table ||), can be better 
understood to consider the Fourier-space version of the Maxwell's equations. 
In fact in many field of physics, including electrodynamics, it is convenient to re-write 
the (Maxwell) differential equations of motion in terms of their Fourier components. In 
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Table I: The transformation properties of the electromagnetic currents, fields and potentials, under 
the discrete symmetries P, T and C, and some important combinations like CP and CPT . If for a 

generic $ of these quantities, including to and k themselves, we write the above transformation rules 

~ x ~ 

in the form $ — > r]x^ , then the value of the multiplicative constants r}x is given for the various 
transformations. The same procedure is followed for the pulsation and wave vector that transforms, for 
example, under the time reversal symmetry, according to (w,k) — > ( w ,k) . 





Transformation 


Quantity 




rjr 


Vc 


Vcp 


rjCPT 




+ 




+ 


+ 




k 




+ 


+ 






P 


+ 


+ 








j 








+ 




E 




+ 




+ 


+ 


B 


+ 








+ 


A 








+ 





this linear approximation, from the usual results of Fourier analysis, we express a general 
solution as a superposition of plane waves through the following relation 

E(i, x) = -—^ J duj J d 3 kE(w,k)exp[i(k-x-^)] (48) 

and similarly for other quantities. Now, if the Fourier transform of a linear equation 
is taken, the space and time drop out. By these means, the set of coupled (Maxwell) 
differential equations are replaced by a set of coupled algebraic equations, which of course 
greatly reduces the difficulty of the problem. The solutions of the algebraic equations are 
the Fourier transforms of the required quantities which can then be obtained by Fourier 
inversion. The Fourier momentum space version of the Maxwell's equations in a medium 
is then given by 

k • B = , kxE = wB (49) 

ik ■ D = p ext , ikxB = fi 3 — i%E (50) 

c z 

with source terms just the external charges and currents. In writing these relations it is 
assumed that the response of the system to the applied electromagnetic fields is linear. 
Even when non linear effects are considered, they could be usually treated as corrections 
to linear ones. In this approach, it is assumed that all quantities which vary with time 
and space differ only slightly from an average value. The source Maxwell's equations (|50| ) 
reduce to the wave form 

^-E + k x (k x E) = -iujuq 3 . (51) 

Eliminating the induced part Ji n d of the current J = Ji n d + J e xt in terms of general 
conductivity tensor a^j 



10 



(jind) . = Vij^j (52) 

then, Eq.(51) reduces to the inhomogeneous wave equation 

A^ = -^(j e:ct )., (53) 

with 



k 2 c 2 ~ 
Aij (w, k) = — r (fcjfcj - <5ij) + tfij (w, k) (54) 

being the equivalent response tensor 

Kij (u, k) = Sij H (Tjj (w, k) . (55) 

If the medium is dissipative, the tensor may be separated into a hermitian and an 
antihermitian part 

Kij = K% + Kf j (56) 

with 

K$ = l(Kij+Kii) ( 57 ) 

which corresponds to a separation into a dissipative and non-dissipative part. The 
Kramers-Kronig relations assure that the hermitian and antihermitian parts are not inde- 
pendent. In fact, 

^ = i,r54i/ (59) 

. r+o° 

KS-Sa = -V -^cLo' . (60) 

For an isotropic medium, that is only spatially dispersive, we have three independent 
(transverse, longitudinal and rotatory) parts in 

e T e R \ 

K$ = e T {5 ij -k i k j ) + e L k i k j +e R e ijr k r = ( -e R e T (61) 

:' / 

where the last matrix form holds if we choose k along the z-axis. The fields E and B 
depend on (u, k), whereas the quantities e T , e L and e R , which represent the properties of 
the medium, are supposed to depend on k only through its magnitude k = |k|. In fact 
we are considering a medium in which only the amplitudes of the electric and magnetic 
polarizations depend on direction. The electromagnetic constants e T , e L and e R have 
no anisotropy and are the same in all the directions. Before we proceed any further, we 
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want to note that there is one relation that e T , e L and e R must satisfy based on very 
general grounds. It is simply the requirement that the fields E, B and the current j are, in 
coordinate space, real quantities. This yields that, for real values of oj, the electromagnetic 
constants satisfy 

r(- w ,k) = /( W ,k), (62) 

where / stands for any of the electromagnetic constants e T , e L and s R and the dependence 
is supposed through its magnitude k, as noted earlier. All that means that, the real parts 
of these electromagnetic constants are even in oj whereas the imaginary parts are odd. We 
stress that this is a very general property, not tied up in any way with the discrete C, P 
and T symmetries of space-time. Of course, for the validity of causality and the consequent 
Kramers-Kronig dispersion relations, oj can in general be allowed to be complex and this 
feature must be included in the analysis. 

The implication of C, P and T symmetries on the electromagnetic constants / imposes 
that the reversed equations result invariant under the C, P and T constraints if 

C invariance =>■ no constraints on e T , e L , and e R , (63) 
P invariance =>■ e R = , (64) 



e T (-oj,k) = e T (oj,k) 
T invariance =>- \ e L (— w,k) = e L (oj,k) , (65) 
e R {-oj,k) = -e R {oj,k) 

and being in general /*(— o>,k) = /(w,k), thus implies that e R must be purely imaginary 
and an odd function of u, whereas e T and e L should be real and even function of oj. In 
the case of a dissipative media, the response tensor components will respect the Onsanger 
relations 

Kij (-oj,k) =Kij (u,k) 
K? j {-u,k)=-K? j {u,k) 
whereas in general (oj, — k) = Kij (oj,k). 

CP invariance obviously gives the same constraint as parity alone does, since charge 
conjugation does not affect. 



CPT invariance 



e T (-oj,k) = e T (oj,k) 

e L {-oj,k) = e L {oj,k) . (67) 
e R (-oj,k) = e R (oj,k) 



From these relations, it is quite clear that the presence of the e^-term implies some 
properties that are asymmetric under P and CP transformations. Similarly, since T and 
CPT imply mutually inconsistent constraints on e R , at least one of them must be violated 
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in order for e to exist, e might be expected to vanish if all interactions conserve C, 
P and T. However since 1957 Q we know that parity is violated by weak interactions 
and since 1964 || we know that CP is also violated although to a much lesser extent. 
On the other hand the induced currents arise because of complicated processes taking 
place within the medium, including (super) weak interactions, so that there is no reason 
why the e^-term should not be present. Of course since the e R -term (if present) violates 
CP it would seem to be extremely tiny. The only constraint on the value of e R seems 
e R (io,k) = e R (—ui,k) which follows from CPT symmetry, since there are strong reasons 
to believe that fundamental interactions between particles can never violate CPT. But 
in the presence of matter, the CPT breaking effects could be relevant due to the effects 
of the medium in the underlying effective theory. So that, CPT asymmetries could be 
large and their invariant constraints become irrelevant. As a result, parity violation in the 
weak interactions of the medium can give rise to a non vanishing value of e R and similarly 
a much smaller T-violating effects should be taken also seriously into account. As we 
found before in Eq. (|65|), if T-invariance holds, e R must be purely imaginary and an odd 
function of u. Then, it is worth noting to the extent that T-violation can be neglected 
in the fundamental interactions, it is conceivable to think a medium where e T and e L 
are constants, (being both even functions of uS), but e R cannot have a non zero constant 
term. In this case a medium becomes dispersive according only to e R . When we consider 
more complicated media, it is possible to obtain P-asymmetries from the medium itself. 
Therefore it is important to consider the passage of electromagnetic radiation through a 
generic medium with a built-in asymmetry. This is the case, for example, of natural sugar 
solution, where one helicity of the radiation is in excess over the other, and physically, it 
manifests itself as the phenomenon of optical activity. 

4 The Optical Analogue of the K° — K° Complex System. 

We seek in electromagnetism an analogue of the two-state mixing problem which leads to a 
non-zero value of an optical parameter which could be associated to the impurity parameter 
e. The problem is quite difficult since the Maxwell's equations are time-reversal invariant. 
The main features of irreversibility enter only considering the effects of anisotropy and 
dissipation. In an isotropic medium, transverse electromagnetic waves, have two degener- 
ate states of polarization, and in anisotropic media, the natural modes are non degenerate 
and not necessarily transverse. Radiation propagating through such a medium splits into 
components in the two natural modes and the difference between the complex refractive 
indices of the two modes produces radiation with interfering effects when the components 
are recombined. In the case of optically active media, the dynamic behavior is connected 
to dissipation more than anisotropy, as far as their electromagnetic properties are con- 
cerned. We review here the fundamental mathematical description of the propagation of 
optical polarization in a medium. In a complex dielectric material, the nature of the wave 
propagation is learned, in general, from the dispersion relations. These are the relations 
that uj and k must satisfy in order to solve the homogeneous wave equation for the E and 
B fields: 

AijEj = . (68) 

In a general complex medium, the existence of a solution is connected to the dispersion 
relation 
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detAy = . (69) 

When this dispersion relation is satisfied, there exists a solution of the homogeneous 
wave equation. In the vacuum, this vector solution E is perpendicular to the direction 
of propagation of the electromagnetic wave. Thus, one can discuss this solution (with 
complex coefficients) 

E = E x e x + E 2 e 2 (70) 

in terms of the components 

E x = E x =Re (E ■ e x ) 

(71) 



E 2 = E y =Re (E ■ e 2 ) 



of the electric field transverse to the coherent beam propagating in the k = z direction. 
In this description of polarization, we refer to two basic polarization orthonormal unit 
vectors e\ and e 2 , both being perpendicular to the wave unit vector k. An equivalent set 
of basic vectors is given by the following right and the left circular basis 

1 / 

f (72) 
e_ =-j= (ei - ie 2 ) . 

Since e 2 = k x ei, these vectors are of the form e\ ± ik x e\. In this language, the effect 
of parity transformation is to reverse the sign of the transverse polarization vector e 2 , 
in a suitable gauge. Whereas, under parity, the circular polarization vectors interchange: 
e_|_ <-> e_. In non dispersive media (e R = 0), the two transverse modes have the same 
physical characteristics, with the same speed of propagation c/ \fejl and, if parity is a good 
symmetry, there would be also no difference in the physical properties of a right-circular 
and left-circular polarization state of the transverse electromagnetic wave. Wherever e R = 
0, in an isotropic homogeneous medium, the response tensor is given by 

Kij = e L k i k i + e T (<% - k^) . (73) 
In this case, the dispersion relation 

detAy = det e L k i k j + (e T - n 2 ) (kjk,- - <%) = e L (e T - n 2 ) 2 = (74) 
yields a longitudinal dispersion relation 

e L {to, k) = (75) 

and a transverse dispersive equation 

n 2 = e 1 (uj,k)= v - = c 2 e (to) (i (to) . (76) 

It is worth noting that only in an isotropic medium, all solutions are polarized strictly 
transverse. However, in an isotropic optically active media (e R / 0), parity, and, in a less 
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extent time-reversal are broken. If also dissipative effects are considered, the two trans- 
verse, (and circular), birefringent (or dichroic) modes have different propagating proper- 
ties. For example, they travel with different speeds through matter and have different 
indices of refraction and coefficients of absorption. It is evident that the two transverse 
modes have different properties when e R does not vanish. The effect of dissipation in 
an isotropic medium will cause a non trivial complex form of e R (u>,k). In fact, if time- 
reversed invariance is satisfied, e T and e L are real, whereas e R is purely imaginary only 
neglecting dissipation. In this case, for a large range of e R , the dispersion relation is solved 
for a real co and k. Thus, the amplitude of the wave will remain unchanged during the 
propagation, which means that there is no absorption in the medium. In the more general 
case of dissipation, time-reversal symmetry can be violated and the space-time evolution of 
the transverse electric field becomes governed by a set of coupled second order differential 
equations 

E 4 = AijEj . (77) 

The two transverse modes no longer have the same dispersion relation. The transverse 
components are governed by the following inverse matrix propagator 

^ = -^r<% + (<% + x*j) ( 78 ) 

which is defined for the transverse E components Eq. (70) with respect to k. However, the 
degree and type of the polarization of this resulting transverse electric wave is associated 



with a unit vector and therefore with a point on the unit sphere [12]. The stereographic 
projection of this unit vector onto the equatorial plane can be consequently associated 
to the ratio of the complex components of a Poincare spinor [|fj|. Nothing else that the 
usual isomorphism between the rotation group and its covering group SU(2). The two- 
components Poincare spinor requires the standard matrix representation in terms of the 
Pauli matrices and it is associated to a vector representing the angular momentum for a 
spin-^ state. In fact, this vector results just the generator of the internal rotations and it is 
called Jones vector [fOfl , when it is considered in the Cartesian basis. Its components give 
the Cartesian components of the complex polarized electric field. This 2x2 unimodular 
matrix representation of the Pauli algebra to describe the polarization, can be replaced by 
the 4x4 Mueller matrix method pq]. Another modern method using 2x2 matrices to 



represent polarization is represented by Fano's density operator approach [16] analogously 
to the case of spin density theory for the spin-i particles. It is worth noting that while the 
fundamental mathematical description of spin I and optical polarization are practically 
the same, their physical interpretations are quite different. In the description of spin i 
we have a spin vector whose direction in a cartesian 3-dimensional space depends upon 
a state vector in a complex 2-dimensional (spin- up, spin-down) spinor space. Anyway, 
in the description of polarization the two ordinary components x and y play the same 
role which the spinor bases played while describing spin ^. Although the problem to 
describe the evolution of the two transverse polarization states of the electromagnetic 
radiation propagating into a complex medium can be described in these several manners, 
nevertheless, a dispersion relation governing the propagation is required. 
In our case of a wave propagating inside a medium with complex e R ^ 0, the space 
propagation of the transverse optical polarization states is governed by the transverse 
response tensor 
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K± = 6 ij + Xij =V 2 (79) 
which is given in terms of a general non Hermitian susceptibility tensor 

x=(T T) ■ ^ 

VX21 Xn J 

Consider a plane polarized wave incident on a medium at z = 0, moving in the positive 
z-direction, the transverse electric components can be assumed in a normal form of the 
type E(z,t) = E{z)e~ wt where 

E(z) = E 1 (z)e 1 + E 2 (z)e 2 (81) 

h-i(S8£) 

where 



i 



= u^JTt fl + xu X12 \ 2 _ ^\fl^ ( 1 + XH + X12 



X21 I + X11/ c(v / AT+v / A2)V X21 1 + X11 + VAiA 2 

(83) 

In the matrix square-root we use the eigenvalues of V 

Ai, 2 = i(E±A) (84) 



with 

E =tr£> 2 = 2 + trx = 2 (1 + X n) 
A 2 =tr 2 £> 2 - 4det£> 2 = tr 2 x - 4detx = 4xi2X2i 
Therefore, the electric field will propagate as 

E + (z) =E + (0)e- lX+z 
E„{z) =E_(0)e- lX ~ z 



(85) 



(86) 



where the V eigenvalues A + and A_ are given by 

A± = A lj2 + 7AlA2~ . (87) 

Here, E + and E- are given by 

E + = /-, ^1 12 ( El + aE ^ 
V 1 + \a\ 



E_ = /T -i 1 -i7 (^1 " ai ^) • 
v 1 + M 

We recognize that the parameter a could be put in relation with the off-diagonal terms of 
the propagation matrix V by the following equation 



>D 21 _ _ K**-Kfi _ 1x21 



A* 



a = 



D l2 V^ii V K* + K& Vxi2 
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It can be viewed that it depends upon the off-diagonal terms of the susceptibility tensor x- 
In a non-dissipative medium K^j = and \a\ = 1. In particular, in an isotropic medium 
which is only dispersive, we have: a = i. The analogy with the neutral kaon system is 
quite evident; the eigenstates the space propagation resemble the mass eigenstates of the 
effective Hamiltonian H governing the evolution of Kl and Ks . 

In summary, it was presented a successful attempt to realize the optical analogue which 
reproduces the peculiar asymmetric effects of the C-P-violating propagation. The ab- 
sorptive and dispersive contributions in isotropic media can then resemble the intriguing 
features of the kaon complex evolution in the realm of high energy particle physics. A 
last remark merits the eventuality to introduce an absorptive effect by means of the space 
anisotropy whose consequent longitudinal propagation, however, seems outside an abelian 
and strictly massless effective theory of the electromagnetic interactions. 
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